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Abstract.— We report patterns of infection of Batrachochytrium dendrobatidis ( Bd) in a cloud forest 
amphibian community in the Andean Western Cordillera of Ecuador (Reserva Las Gralarias). Data 
were obtained during the rainy seasons of two consecutive years, using qPCR (year 2012) and end¬ 
point PCR (year 2013). We show that average Bd prevalence in this amphibian community is high 
(2012:35-49%; 2013:14-32%), but found no evidence of population declines or that Bd is negatively 
affecting host populations. We found a significant correlation between Bd prevalence and taxonomy, 
reproductive mode, and habitat, but no correlation between Bd infection intensity and the same three 
variables. Contrary to our expectations, frog species with aquatic reproductive modes (glassfrogs, 
treefrogs) showed lower Bd prevalence than direct-developing frogs ( Pristimantis spp.). Although 
further monitoring is needed to determine long-term population trends, our two-year dataset on 
disease and population size support the hypothesis that frogs are tolerant to infection, a condition 
that could potentially have resulted from exposure to previous Bd epidemic outbreaks. 

Resumen. —En este estudio reportamos datos sobre los patrones de infeccion de Batrachochytrium 
dendrobatidis ( Bd) en una comunidad de anfibios en la Cordillera Occidental de los Andes del 
Ecuador (Reserva Las Gralarias). Los datos fueron obtenidos durante la estacion lluviosa en dos 
anos consecutivos, utilizando qPCR (ano 2012) y PCR de punto final (afio 2013). Los resultados 
muestran una alta prevalencia de Bden la comunidad (2012:35-49%; 2013:14-32%); sin embargo, no 
se encontro evidencia de disminuciones poblacionales o de que Bd este afectando negativamente a 
las especies de anfibios. Existe una relacion significativa entre la prevalencia de Bdy la taxonorma, 
modo reproductive y habitat de los anfibios, pero no hubo correlacion entre la intensidad de 
infeccion de Bd y las mismas tres variables. Contrario a nuestras predicciones, las especies de 
anuros con larvas acuaticas (ranas de cristal, ranas arboreas) presentaron prevalencias de Bd mas 
bajas que los anuros de desarrollo directo ( Pristimantis spp.). A pesar de que se requiere de un 
monitoreo continuo para determinar las dinamicas poblacionales a largo plazo, los datos obtenidos 
hasta el momento apoyan un escenario donde las especies de ranas de la Reserva Las Gralarias 
parecen tolerar la infeccion de Bd, una condicion posiblemente adquirida mediante la exposicion a 
brotes epidemicos previos. 
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Introduction 

A third of global amphibian species are threatened with 
extinction (Stuart et al. 2004; Wake and Vredenburg 
2008) and, most concerning, numerous local popula¬ 
tion declines and extinctions have occurred in relatively 
pristine areas, where anthropogenic habitat destruction is 
low (Lips 1998, 1999; Drost and Fellers 1996; La Marca 
et al. 2005). In the last two decades, several studies have 
attributed anuran mass mortality events to the emergence 
of the pathogenic fungus Batrachochytrium dendroba- 
tidis (Bd ), a pathogen with widespread geographic and 
ecological distribution (Berger et al. 1998; Daszak et al. 
1999, 2003; Lips et al. 2006; Becker and Zamudio 2011; 
Rodriguez et al. 2014). Alternative explanations to am¬ 
phibian declines add a role to global warming and tem¬ 
perature variability (Pounds et al. 2006; Rohr and Raffel 
2010; Menendez-Guerrero and Graham 2013). 

Batrachochytrium dendrobatidis infects the keratin¬ 
ized skin of amphibians and disrupts the regulatory func¬ 
tioning of the integument (Berger et al. 1998; Voyles et 
al. 2009). Infection inhibits host immune responses in 
some species (Fites et al. 2013) and in severe cases of 
infection, electrolyte depletion and osmotic imbalance 
may lead to mortality (Voyles et al. 2007, 2009). How¬ 
ever, not all amphibian species are equally susceptible to 
the pathogen. For example, at Santa Fe, Panama, Bd has 
caused declines or local extinctions of most anurans in 
the original community, but six species of frogs and toads 
remain abundant, despite being infected by the fungus 
(Lips et al. 2006). In laboratory challenge experiments, 
amphibian mortality rates range from 0% to 100%, de¬ 
pending on the species, host age, pathogen genotype, and 
dosage (Berger et al. 2005a; Daszak et al. 2004; Longo et 
al. 2014). The reasons for host differences in susceptibil¬ 
ity include immunogenic variation (Ellison et al. 2014; 
Savage et al. 2014), microhabitat use (Kriger and Hero 
2007; Grundler et al. 2012), association with water as 
embryos, tadpoles, or adults (Lips et al. 2003), and host 
thermoregulatory behavior (Richards-Zawacki 2010). 
Because Bd transmission may happen through frog-frog 
contact, or through motile zoospore movement from one 
host to another, frogs and toads that spend more time in 
water are expected to have higher exposure and suscep¬ 
tibility to infection than species that are primarily terres¬ 
trial (i.e, direct developers; Lips et al. 2003; Kriger and 
Hero 2007). 

In this study, we report data on infection patterns of 
Bd obtained during the rainy seasons of two consecu¬ 
tive years in the amphibian community of Reserva Las 
Gralarias, a cloud forest site in the Andean Western Cor¬ 
dillera of Ecuador. We found that Bd prevalence in all 
amphibian species is high, but found no evidence that 
Bd is negatively affecting amphibians (i.e., no apparent 
population declines, or records of clinical signs of chy- 
tridiomycosis). We report on infection prevalence and 
intensity dynamics for the two-year period, and examine 
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patterns of Bd infection in species that vary in their tax¬ 
onomy, reproductive mode, and habitat. 

Materials and Methods 

Study site: The study was conducted at Reserva Las 
Gralarias (0°01’ S, 78°44’ W; 1822-2400 m), a private 
reserve covering an area of 1,063 acres (425 ha) located 
on the Pacific slopes of the Andes, Pichincha Province, 
Mindo Parish, Ecuador. The study site has an eleva- 
tional range of 1,825-2,400 m and includes primary and 
secondary forest, regenerating pasture, and numerous 
ephemeral and permanent streams and creeks (Hutter and 
Guayasamin 2012). 

Amphibian taxonomy: For generic and suprageneric clas¬ 
sification, we follow the taxonomic proposals of Hedges 
et al. (2008), Guayasamin et al. (2009), and Faivovich et 
al. (2005), Pyron and Wiens (2011), as summarized in 
Frost (2014). 

Amphibian richness and abundance: During the rainy 
seasons of 2012 (23 January-29 March) and 2013 (14 
March-22 April), we sampled trails of Reserva Las 
Gralarias during the night, including most of its habitat 
heterogeneity and elevational gradient, to record the spe¬ 
cies richness of the reserve. We placed eight transects, 
each with an area of 500 x 4 m (Appendix 1), to maxi¬ 
mize species detection and to obtain a baseline dataset 
on population size and Bd prevalence. Each transect was 
sampled by two people for 3-4 hours during the night 
(generally starting at 8 pm); temperatures during sam¬ 
pling varied between 11-15 °C. All detected amphibians 
were, when possible, photographed. Calling males were 
also reported and identified with the aid of photographic 
and acoustic guides (Arteaga et al. 2013; Centro Jambatu 
2011-2014). We used a Student’s t-test to quantify differ¬ 
ences in population sizes in transects that were sampled 
multiple times during the rainy season of 2012 and 2013 
(Lucy’s Creek and Kathy’s Creek); the normality of spe¬ 
cies abundance was assessed using a Shapiro-Wilk Test. 

Diagnosis of Batrachochytrium dendrobatidis: We 
swabbed the ventral regions of all amphibians captured 
in our survey, following the standard procedures in Hyatt 
et al. (2007; Fig. 1); dry swabs were stored in -4 °C until 
analysis. Testing for Bd was carried out using Real-Time 
PCR (q-PCR) for samples obtained during 2012 and end¬ 
point Polymerase Chain Reaction (PCR) for samples ob¬ 
tained in 2013; the use of these two methods was contin¬ 
gent on access to q-PCR (available during 2012). In both 
cases, DNA extractions were carried out using guanidin- 
ium thiocyanate. For samples obtained during 2012, we 
used a 1:10 dilution of the extract as template in Taqman 
q-PCR assays for the detection of Bd (Boyle et al. 2004). 
This assay uses 5d-specific primers ITS 1-3 Chytr and 
5.8S Chytr, in addition to the fluorescently-labeled probe 
Chytr MGB2, and amplifies the ITS-1 fragment of the 

August 2014 | Volume 8 | Number 1 | e81 


Batrachochytrium dendrobatidis in an Andean frog community 



Fig. 1. Swab sample obtained from Centrolene heloderma at 
Reserva Las Gralarias, Ecuador. 


Bd genome at the junction of the ITS-1 and 5.8S regions. 
We used a standard curve that included 1000, 100, 10, 
1, and 0.1 zoospore genome equivalents, and followed 
qPCR conditions described in Boyle et al. (2004). For 
samples obtained during 2013, Bd presence was tested 
using the internal transcribed spacer regions (ITS-1, ITS- 
2) primers Bdla (5*-CAGTGTGCCATATGTCACG-3’) 
and Bd2a (5 * -CATGGTTCATATCTGTCCAG-3 ’) de¬ 
veloped by Annis et al. (2004); the presence/absence of 
Bd was determined via the visualization of the amplified 
band in agarose gel electrophoresis. The two methods to 
detect Bd have different sensitivities; therefore, direct 
comparisons of Bd prevalence between years should be 
considered with caution. However, family and habitat 
correlates with infection status should not be biased by 
detection method, and qPCR offers the additional advan¬ 
tage of quantifying infection intensity (load). 

Prevalence and correlates of Batrachochytrium den¬ 
drobatidis in amphibians: We estimated prevalence of 


Bd within each anuran species as the number of frogs 
that tested positive for Bd, divided by the total number of 
sampled frogs for that particular species in a given year. 
We estimated the 95% confidence interval for preva¬ 
lence in each species, (Wilson 1927; Newcomb 1998). 
We modeled Bd presence or absence in each individual 
by using a logistic regression. We tested for possible as¬ 
sociations of Bd prevalence with the following variables: 
habitat (terrestrial, riparian, lentic), reproductive mode 
(aquatic, terrestrial), and taxonomy (family). Statistical 
significance of results was assessed with a chi-square 
test. 

Infection intensity of Batrachochytrium dendrobati¬ 
dis and correlates in amphibians: We tested for possible 
associations of Bd intensity (measured as zoospore ge¬ 
nomic equivalents) with the following variables: habitat 
(terrestrial, riparian, lentic), reproductive mode (aquatic, 
terrestrial), and taxonomy (family, genus, species; Ap¬ 
pendix 2). Given the strong right skew of infection load, 
we used the non-parametric Kruskal-Wallis test. All sta¬ 
tistical analyses were performed using R v. 2.15.3 (R 
CoreTeam 2012). 

Results 

Species richness and abundance: During the two sam¬ 
pling periods, we recorded a total of 2,450 individuals 
of 28 species (Appendix 2). The abundance of species at 
Lucy’s Creek and Kathy’s Creek is summarized in Tables 
1 and 2. Because most taxa were scarce, we restricted 
the comparisons between years to relatively abundant 
species (glassfrogs). Abundances of glassfrogs at Lucy’s 
Creek and Kathy’s Creek were not significantly different 
between years (Tables 1, 2). 

Prevalence of Batrachochytrium dendrobatidis in 
amphibians: Swabs of 320 frogs were tested for Bd, and 


O 



Fig. 2. Significant Bd infection differences in amphibians according to reproductive modes, habitat use, and taxonomy. P values are 
reported for 2012 and 2013; significance is noted by ** (p < 0.01) and *** (p < 0.001). 
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Table 1 . Abundance of amphibian species at Lucy’s Creek, Reserva Las Gralarias. Abundances are presented as minimum-maxi¬ 
mum, followed, in parenthesis, by mean + standard error. The Student’s /-test was performed only in species with normally distrib¬ 
uted abundances. 



Lucy’s Creek 

f-test (p) 

Population trend 

Year 

2012 

2013 



Number of nights sampled 

10 

6 



Family/Species 

Centrolenidae 





Centrolene lynchi 

5-20(13.1+5.13) 

2-18 (9.0 + 6.54) 

0.184 

No difference 

Centrolene peristictum 

4-35 (20.7 + 11.68) 

6-26(15.2 + 7.00) 

0.314 

No difference 

Nymphargus grandisonae 

1-10 (6.2 ±3.12) 

0-8 (3.8 + 2.93) 

0.155 

No difference 

Hylidae 





Hyloscirtus alytolylax 

3-7 (4.3 + 1.42) 

0-6(3.3 + 2.16) 

0.295 

No difference 

Craugastoridae 





Pristimantis achatinus 

0-2 (0.3 + 0.67) 

0-1 (0.2 + 0.41) 

— 

— 

Pristimantis appendiculatus 

0-2(0.8+0.92) 

0-1 (0.5 + 0.55) 

— 

— 

Pristimantis calcarulatus 

0-2 (0.2 + 0.63) 

0-1 (0.2 + 0.41) 

— 

— 

Pristimantis eremitus 

0 

0-2 (0.5 + 0.84) 

— 

— 

Pristimantis eugeniae 

0 

0-2 (0.3 + 0.82) 

— 

— 

Pristimantis illotus 

0-1 (0.1+0.32) 

0 

— 

— 

Pristimantis parvillus 

0-1 (0.1 +0.32) 

0 

— 

— 

Pristimantis sobetes 

0 

0-2 (0.3+0.82) 

— 

— 

Pristimantis w-nigrum 

0-2 (0.6 + 0.84) 

0-2(0.7 + 0.82) 

— 

— 


approximately a third of those were positive. In samples 
from 2012, prevalence of Bd was relatively high, with 
42% of all frogs testing positive for Bd infection. Dur¬ 
ing 2013, Bd prevalence was 22%. Differences in preva¬ 
lence between the two years are likely caused by detec¬ 
tion method. Most species infected in 2012 carried low 
Bd loads as determined by qPCR; the highest Bd load 
obtained was in Centrolene ballux with 22.5 genomic 
equivalents. Prevalence per species per year is summa¬ 
rized in Table 3. 

The logistic regression shows a significant relation¬ 
ship (p < 0.001) of Bd infection with species reproduc¬ 
tive mode, habitat, and taxonomy (Fig. 2). Frogs with 
a terrestrial reproductive mode (direct developers; i.e., 
genus Pristimantis ; see Duellman and Trueb 1986) have 
a higher Bd prevalence than amphibians with aquatic re¬ 
production (i.e., glassfrogs and treefrogs). Frog species 
that are dependent on riverine habitats for reproduction 
show significantly less infection than anurans that use 
terrestrial or lentic habitats for reproduction (p < 0.001). 
Also, species in the Centrolenidae family (glassfrogs) 
show a lower Bd prevalence than species in Craugastori- 
dae and Hylidae (Table 3). Although, Bd prevalence dur¬ 
ing 2012 was significantly higher than in 2013 (probably 
as a result of higher sensitivity of qPCR), we found no 
significant interaction among sampling year and repro¬ 
ductive mode, habitat, or taxonomy. 

Infection intensity of Batrachochytrium dendrobati- 
dis and correlates in amphibians: We found no relation¬ 


ship between Bd infection intensity (Bd load, year 2012) 
and taxonomy, reproductive mode, or habitat. 

Discussion 

Our results show a relatively high mean prevalence of 
Bd (36%) across both years in the Andean frog commu¬ 
nity of Reserva Las Gralarias (see Hossack et al. 2010 for 
comparison). From a total of 20 species analyzed, only 
three (Nymphargus griffithsi , Pristimantis illotus , and P 
pteridophilus ) tested negative for Bd; however, sample 
sizes for non-infected species were low (5, 1, and 3 in¬ 
dividuals, respectively). Infected species included frogs 
with very different reproductive modes, including taxa 
with terrestrial direct development (Pristimantis spp.), 
species that deposit eggs in ponds (Dendropsophus car- 
nifex ), and others that place their eggs on vegetation from 
where hatching tadpoles drop into streams (Centrolene 
spp., Nymphargus spp., Hyloscirtus spp.). 

Because Bd is an aquatic pathogen (Berger et al. 
2005b) we expected amphibian species with aquatic re¬ 
productive modes to show higher infection prevalence 
(Lips et al. 2005; Brem and Lips 2008). In fact, the most 
dramatic amphibian declines and extinctions in the Andes 
have occurred in species with aquatic larvae (La Marca 
et al. 2005; Bustamante et al. 2005; Merino-Viteri et al. 
2005; Coloma et al. 2010). Our results indicate, surpris¬ 
ingly, a higher Bd prevalence in frogs with a terrestrial 
reproductive mode (Pristimantis spp.) than in those that 
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Table 2. Abundance of amphibian species at Kathy’s Creek, Reserva Las Gralarias. Abundances are presented as minimum-maxi¬ 
mum, followed, in parenthesis, by mean + standard error. The Student’s t -test was performed only in species with normally distrib¬ 
uted abundances. 



Kathy’s Creek 

f-test 

Population trend 

Year 

2012 

2013 



Number of nights sampled 

10 

5 



Family/Species 

Centrolenidae 





Centrolene ballux 

5-37 (22.7+ 11.6) 

3-25 (11.4 + 8.67) 

0.078 

No difference 

Centrolene peristictum 

0-5 (2.1 + 1.66) 

0-5 (2.0+1.87) 

0.918 

No difference 

Nymphargus grandisonae 

0-7(3.7 + 2.21) 

0-6 (2.2 + 2.28) 

0.242 

No difference 

Nymphargus griffithsi 

0-8 (2.3 + 2.26) 

0-3 (1.4+1.34) 

— 

— 

Nymphargus lasgralarias 

3-28 (19.4 + 8.54) 

7-28 (15.0 + 8.69) 

0.366 

No difference 

Hylidae 

Hyloscirtus alytolylax 

Craugastoridae 

0-1 (0.1 +0.32) 

0-1 (0.4 + 0.59) 

— 

— 

Pristimantis achatinus 

0-3 (0.3 + 0.95) 

0-1 (0.2 + 0.45) 

— 

— 

Pristimantis appendiculatus 

0-7(1.2 + 2.10) 

0-1 (0.4 + 0.59) 

— 

— 

Pristimantis calcarulatus 

0-3(1.1+0.74) 

0-3 (1.0+ 1.23) 

— 

— 

Pristimantis eremitus 

0-1 (0.1+0.32) 

0-1 (0.2 + 0.45) 

— 

— 

Pristimantis eugeniae 

0-1 (0.1 +0.32) 

1-2 (0.8 + 0.84) 

— 

— 

Pristimantis sobetes 

0-1 (0.2 + 0.42) 

0-2 (0.4 + 0.89) 

— 

— 

Pristimantis w-nigrum 

0-1 (0.1 +0.32) 

0-1 (0.4 + 0.59) 

— 

— 


reproduce in water (mainly glassfrogs; Centrolene spp., 
Nymphcirgus spp.). This finding supports the idea that 
even terrestrial breeders may serve as reservoirs for the 
pathogen in diverse amphibian communities (Longo et 
al. 2013). Higher prevalence in terrestrial frogs requires 
that Bd zoospores survive in terrestrial habitats. John¬ 
son and Speare (2003) indicated that Bd can survive in 
moist soil for up to three months. Cloud forests in west¬ 
ern Ecuador typically have near constant rain and high 
levels of humidity during the rainy season (Hutter and 
Guayasamin 2012; Arteaga et al. 2013), and this may ex¬ 
tend zoospore survival in terrestrial environments at Las 
Gralarias. Higher Bd prevalence in terrestrial frogs com¬ 
pared to that in frogs with aquatic reproduction might 
also be related to intrinsic differences in, for example, the 
efficacy of immune responses (Rosenblum et al. 2009; 
Woodhams et al. 2007) or differences in anuran skin mi¬ 
crobiota (Flechas et al. 2012). 

Our results also show that Bd prevalence is signifi¬ 
cantly associated with taxonomy (i.e., family). Thus, 
glassfrogs (family Centrolenidae) might have immune 
responses or skin microbiota that work as better barri¬ 
ers to the pathogen than those in terrestrial (i.e., Pristi- 
mcintis ) frogs. The strong correlation of prevalence with 
taxonomy, habitat, and reproductive mode (Appendix 2) 
indicates that further studies need to focus on the specific 
effects of each of these factors; in other words, phylog- 
eny (and taxonomy) correlates with reproductive mode 
and habitat use. 

A second surprising finding of our study is that, al¬ 
though prevalence of Bd is high in most anuran species, 


we did not observe any sign of population declines or 
abrupt crashes, nor have we found dead or sick frogs dur¬ 
ing four years of intensive fieldwork (2010-2014; JMG 
pers. obs.). Thus, this frog community persists with an 
endemic pathogen and with relatively low loads (less 
than 10 zoospores, Table 3). The apparent increased re¬ 
sistance or tolerance of amphibians from Reserva Las 
Gralarias to Bd infection may be explained by one or 
several of the following mechanisms: (i) amphibian in¬ 
nate and/or acquired defense mechanisms (Savage and 
Zamudio 2012; Woodhams et al. 2007); (ii) skin bacte¬ 
rial commensals with anti-fungal properties (Harris et al. 
2006); (iii) behavioral and ecological factors that reduce 
the likelihood of infection and disease (e.g., microhabitat 
selection, reproductive mode; Lips et al. 2003; Rowley 
and Alford, 2007), and/or (iv) variation in Bd virulence 
(Berger et al. 2005a; Fisher et al. 2009). Our amphib¬ 
ian monitoring took place during the rainy season when 
most species are active. Therefore, future studies should 
determine whether or not this apparent tolerance to Bd 
is stable through longer periods of time or if it fluctuates 
depending on environmental variables influencing host 
immunity, behavior, microbiota, or pathogenicity. 

After the emergence of an infectious disease, surviv¬ 
ing hosts can evolve tolerance or resistance (Retallick et 
al. 2004; Savage and Zamudio 2011). The earliest known 
record of Bd in Ecuador is in 1980, in the Harlequin frog 
Atelopus bomolochos (Ron and Merino-Viteri 2000), 
a species that is now probably extinct (Coloma et al. 
2014). If Bd reached and spread in Ecuador during the 
early 1980s (Ron et al. 2003; Lips et al. 2008), we hy- 
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Table 3. Prevalence of Batrachochytrium dendrobatidis ( Bd ) in amphibians at Reserva Las Gralarias, Ecuador, during the rainy 
seasons of 2012 and 2013. Bd prevalence for each species is followed, in parenthesis, by a 95% confidence interval. Bd load sum¬ 
mary data includes only samples that tested positive for Bd. 



2012 

2013 


N 

Positive Bd 

Prevalence 

Bd Load 
(mean ± sd) 

N 

Positive Bd 

Prevalence 

Family: Centrolenidae 








Centrolene ballux 

17 

8 

47% 

(24-71%) 

6.5 ± 10.7 

9 

2 

22% 

(39-59%) 

Centrolene heloderma 

6 

1 

17% 

(1-63%) 

0.6 

1 

0 

0% 

(0-95%) 

Centrolene lynchi 

6 

1 

17% 

(1-63%) 

- 

5 

1 

20% 

(1-70%) 

Centrolene peristictum 

21 

6 

29% 

(12-52%) 

2.1 ± 1.5 

16 

3 

19% 

(5-46%) 

Nymphargus grandisonae 

21 

5 

24% 

(9-48%) 

4.7 

- 

- 

- 

Nymphargus griffithsi 

3 

0 

0% 

(1-69%) 

- 

2 

0 

0% 

(0-80%) 

Nymphargus lasgralarias 

16 

6 

38% 

(16-64%) 

6.5 ± 1.5 

10 

2 

20% 

(4-56%) 

Family: Hylidae 








Dendropsophus carnifex 

10 

5 

50% 

(20-80%) 

- 

- 

- 

- 

Hyloscirtus alytolylax 

9 

8 

89% 

(51-99%) 

2.4 ± 2.4 

7 

2 

29% 

(5-70%) 

Family: Craugastoridae 








Pristimantis achatinus 

7 

4 

57% 

(20-88%) 

- 

- 

- 

- 

Pristimantis appendiculatus 

23 

10 

44% 

(24-65%) 

1.9 ±2.4 

- 

- 

- 

Pristimantis calcarulatus 

15 

2 

13% 

(2-42%) 

1.1 

15 

1 

7% 

(4_34%) 

Pristimantis eremitus 

4 

4 

100% 

(40-100%) 

0.9 ±0.2 

4 

2 

50% 

(9-91%) 

Pristiman tis eugeniae 

18 

12 

66% 

(41-86%) 

2.5 

2 

1 

50% 

(3-97%) 

Pristimantis hectus 

8 

2 

25% 

(4-64%) 

- 

14 

4 

29% 

(10-58%) 

Pristimantis illotus 

- 

- 

- 

- 

1 

0 

0% 

(0-95%) 

Pristimantis parvillus 

9 

4 

44% 

(15-77%) 

- 

- 

- 

- 

Pristimantis sobetes 

8 

3 

38% 

(10-74%) 

- 

9 

3 

33% 

(9-69%) 

Pristimantis pteridophilus 

- 

- 

- 

- 

3 

0 

0% 

(0-69%) 

Pristimantis w-nigrum 

21 

13 

62% 

(39-81%) 

1.5 ±0.7 

- 

- 

- 

TOTAL 

222 

94 

42% 

(35-49%) 


88 

19 

22% 

(14-32%) 


pothesize that many of the population declines observed 
in the country at that time (e.g., Coloma 1995, 2002; 
Coloma et al. 2000; Ron et al. 2003; Bustamante et al. 
2005; La Marca et al. 2005; Lips et al. 2008; Coloma et 
al. 2010) could be attributable to chytridiomycosis. Thus, 
it is probable that most Andean amphibian communities 


have been exposed to Bd for more than three decades and 
that current sampling finds remnant species that are tol¬ 
erant to Bd while the susceptible species are already ex¬ 
tinct. Under this scenario, selection should have favored 
the persistence of amphibian species or specific popula¬ 
tions that have developed defenses against Bd; therefore, 
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changing host composition of these communities right 
after pathogen emergence. We acknowledge, however, 
that this is a working hypothesis which assumptions 
depend on our knowledge of the historical distribution 
of the chytrid. For example, if new data shows that Bd 
was present in the Ecuadorian Andes before amphibian 
declines were noticed, such piece of information would 
support the endemic pathogen hypothesis, which states 
that environmental changes triggered Bd outbreaks (Ra- 
chowicz et al. 2005). 

Reserva Las Gralarias is one of the most studied and 
species-rich area in the cloud forest of the tropical An¬ 
des, containing numerous species considered endangered 
by the IUCN (2014; see Appendix 2). However, the 
community (and surrounding areas) lacks at least three 
groups of species that were conspicuous in Ecuadorian 
cloud forests: marsupial frogs ( Gcistrotheca plumbea , 

G. guentheri ), harlequin frogs ( Atelopus longirostris , A. 
mindoensis ), and dendrobatid frogs ( Hyloxalus lehmani , 

H. maquipucuna ) (Coloma et al. 2011-2014; Arteaga et 
et al. al. 2013). Marsupial and harlequin frogs are par¬ 
ticularly susceptible to Bd (Lips et al. 2003; Flechas et al. 
2012; Ellison et al. 2014; DiRenzo et al. 2014) and are 
the primary species that suffered population declines and 
extinctions in Ecuador (Lips et al. 2002; La Marca 2005) 
even in pristine areas. The absence of these lineages at 
Reserva Las Gralarias supports to the hypothesis that this 
is a post-decline amphibian community. Understanding 
the long-term effects of pathogens (eg., chytrid) and tem¬ 
perature variability in such a community is essential for 
the continued effective management of endangered spe¬ 
cies in the Andean cloud forests. 

Considering Ecuador’s high diversity of amphibian 
species, life history modes, and evolutionary history, our 
study provides a baseline to study the evolution of de¬ 
fense strategies against Bd. We recommend further re¬ 
search to determine the mechanisms driving the observed 
differences in pathogen exposure among hosts differing 
in reproductive modes, habitat, and taxonomy. 
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Appendix 1. Transects sampled at Reserva Las Gralarias. Each transect has an area of 500 x 4 m. Latitude and longitude are in 
decimal degrees. 


Transect 

Elevation (m) 

Latitude 

Longitude 

Habitat 

Lucy’s creek 

1822-1858 

At start: -0.00492 

At end: -0.00342 

At start: -78.73344 

At end: -78.74051 

Riverine vegetation along creek 

Kathy’s creek 

2041-2066 

At start: -0.01696 

At end: -0.0156 

At start: -78.7314 

At end: -78.73386 

Riverine vegetation along creek 

Santa Rosa river 

1884-1882 

At start: -0.0133 

At end: -0.01054 

At start: -78.72368 

At end: -78.7211 

Riverine vegetation along river 

Waterfall trail 

1897-2107 

At start: -0.0135 

At end: -0.01379 

At start: -78.72461 

At end: -78.7269 

Primary and secondary terra firme forest 

Five-Frog creek 

2141-2156 

At start: -0.03166 

At end: -0.03098 

At start: -78.70421 

At end: -78.70853 

Riverine vegetation along creek 

Osoverde & Guarumo trail 

2141-2156 

At start: -0.03166 

At end: -0.03098 

At start: -78.70421 

At end: -78.70853 

Primary and secondary terra firme forest 

Puma trail 

1923-2031 

At start: -0.00954 

At end: -0.00708 

At start: -78.7346 

At end: -78.73662 

Primary and secondary terra firme forest 

TKA trail 

2192-2216 

At start: -0.0275 

At end: -0.02516 

At start: -78.70477 

At end: -78.70353 

Primary and secondary terra firme forest 

Peccary trail 

1803-1896 

At start: -0.00750 

At end: -0.0076 

At start: -78.72635 

At end: -78.72862 

Primary and secondary terra firme forest 


Appendix 2. Amphibians at Reserva Las Gralarias, with corresponding IUCN (2014) conservation status. The list includes three 
potential new species ( Pristimantis sp. 1, Pristimantis sp. 2, and Pristimantis sp. 3). Reproductive modes are sensu Haddad and 
Prado (2005). 


Species 

Reproductive mode 

Habitat for reproduction 

Conservation status 

Family: Centrolenidae (7 spp.) 

Mode 25: Eggs hatching into 
exotrophic tadpoles that drop in 
lotic water 



Centrolene ballux 

Mode 25 

Vegetation along fast-flowing 
streams 

Critically Endangered 

Centrolene heloderma 

Mode 25 

Vegetation along fast-flowing 
streams 

Critically Endangered 

Centrolene lynchi 

Mode 25 

Vegetation along fast-flowing 
streams 

Endangered 

Centrolene peristictum 

Mode 25 

Vegetation along fast-flowing 
streams 

Vulnerable 

Nymphargus griffithsi 

Mode 25 

Vegetation along fast-flowing 
streams 

Vulnerable 

Nymphargus grandisonae 

Mode 25 

Vegetation along fast-flowing 
streams 

Least Concern 

Nymphargus lasgralarias 

Mode 25 

Vegetation along fast-flowing 
streams 

Data Deficient 
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Appendix 2 (continued). Amphibians at Reserva Las Gralarias, with corresponding IUCN (2014) conservation status. The list 
includes three potential new species (Pristimantis sp. 1, Pristimantis sp. 2, and Pristimantis sp. 3). Reproductive modes are sensu 
Haddad and Prado (2005). 


Species 

Reproductive mode 

Habitat for reproduction 

Conservation status 

Family: Craugastoridae (16 
spp.) 

Mode 23: Direct development of 
terrestrial eggs 



Pristimantis achatinus 

Mode 23 

Terrestrial, mainly in pastures and 
modified environments 

Least Concern 

Pristimantis appendiculatus 

Mode 23 

Terrestrial, mainly primary and 
secondary forests 

Least Concern 

Pristimantis calcarulatus 

Mode 23 

Terrestrial, mainly primary and 
secondary forests 

Vulnerable 

Pristimantis crenunguis 

Mode 23 

Terrestrial, mainly primary and 
secondary forests 

Endangered 

Pristimantis eremitus 

Mode 23 

Terrestrial, mainly primary and 
secondary forests 

Vulnerable 

Pristimantis eugeniae 

Mode 23 

Terrestrial, mainly primary and 
secondary forests 

Endangered 

Pristimantis liectus 

Mode 23 

Terrestrial, mainly primary and 
secondary forests 

Data deficient 

Pristimantis illotus 

Mode 23 

Terrestrial, mainly primary and 
secondary forests 

Near Threatened 

Pristimantis pan’illus 

Mode 23 

Terrestrial, mainly primary and 
secondary forests 

Least Concern 

Pristimantis pteridophilus 

Mode 23 

Terrestrial, mainly primary and 
secondary forests 

Endangered 

Pristimantis sobetes 

Mode 23 

Terrestrial, mainly primary and 
secondary forests 

Endangered 

Pristimantis verecundus 

Mode 23 

Terrestrial, mainly primary and 
secondary forests 

Vulnerable 

Pristimantis w-nigrum 

Mode 23 

Terrestrial, mainly primary and 
secondary forests 

Least Concern 

Pristimantis sp. 1 

Mode 23 

Terrestrial, mainly primary and 
secondary forests 

Not evaluated 

Pristimantis sp. 2 

Mode 23 

Terrestrial, mainly primary and 
secondary forests 

Not evaluated 

Pristimantis sp. 3 

Mode 23 

Terrestrial, mainly primary and 
secondary forests 

Not evaluated 

Family: Hylidae (3 spp) 

Dendropsophus carnifex 

Mode 1: Eggs and exotrophic 
tadpoles in lentic water 

Ponds 

Least Concern 

Hyloscirtus alytolylax 

Mode 25 

Vegetation along fast-flowing 
streams 

Near Threatened 

Hyloscirtus criptico 

Mode 25 

Vegetation along fast-flowing 
streams 

Not evaluated 

Family: Caeciliidae (1 sp.) 

Caecilia buckleyi 

Family: Rhinatrematidae (1 sp.) 

Unknown 

Unknown 

Not evaluated 

Epicrionops bicolor 

Unknown 

Unknown 

Least Concern 
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